Finding Traps in Non-linear Spin Arrays 
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Precise knowledge of the Hamiltonian of a system is a key to many of its applications. Tasks 
such state transfer or quantum computation have been well studied with a linear chain, but hardly 
with systems, which do not possess a linear structure. While this difference does not disturb the 
end-to-end dynamics of a single excitation, the evolution is significantly changed in other subspaces. 
Here we quantify the difference between a linear chain and a pseudo-chain, which have more than 
one spin at some site (block). We show how to estimate a number of all spins in the system and the 
intra-block coupling constants. We also suggest how it is possible to eliminate excitations trapped 
in such blocks, which may disturb the state transfer. Importantly, one uses only at-ends data and 
needs to be able to put the system to either the maximally magnetized or the maximally mixed 
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state. This can obtained by controlling a global decoherence parameter, such as temperature. 


A very interesting application of spin chains, and more 
generally, lattices, is the quantum state transfer. Bose 
noticed that having a fully magnetized state in the þe- 
ginning, a spin flip propagates from the sender’s site all 
over the system. One may then use a protocol which does 
not require periodic dynamics, but relies on extended in- 
frastructure [2], or a time-dependent unitary transforma- 
tion ia to retrieve the whole information. Alternatively, 
one could use a system with a specific periodic evolution, 
which mirrors the information with respect to the mid- 
dle of the chain. Such a transfer was first introduced in 
id, and independently in B, lel, and the state mirroring 
condition was formalized by Shi et al. [7]. 


In particular, some research has been devoted to the 
state transfer in the regime of limited access to the mid- 
dle of the chain. This is a natural assumption, as in 
the microscale the quantum link between users is a black 
box. It was shown, for example, that with an rx spin-4 
a quantum state can be transferred without initializing 
the state of the interconnecting part of the chain |8], and 
if users can manipulate two spins, it can be even done 
without remote collaboration (gj (see also id). Inter- 
estingly, a toy model analysis suggests that this feature 
allows to better transfer information in presence of some 
types of decoherence {1}. It was also shown that limited 
access still allows to perform a Hamiltonian tomography 
of whether xx or zzz chains and more compli- 
cated lattices |14]. It has also been shown that acting 
on extreme sites of a chain one can perform quantum 
computation [16]. 


In this contribution we present a simple toolbox, 
which allows to estimate tomography of a pseudo-chain. 
The end-to-end evolution in the one excitation subspace 
(where only one spin is oriented up rather than down, 
OES) governed by a pseudo-chain is exactly the same 
as of the corresponding linear chain, but it is differ- 
ent in other subspaces. A pseudo-chain has a block 
of spins instead of a single one at some site. These 


spins have all equal coupling strengths to neighboring 
sites (blocks). Considering OES evolution, if these cou- 
plings were J;, the effective coupling between blocks is 
VJ NiNi+idi, Ni being the number of physical spins in 
block i (Ni = Nw = 1). An illustrative example is 
a map between Christandl et al.-Nikopoulos-Petrosyan- 
Lambropoulos chain and a hypercube of equally coupled 
spins described in ah. This might be of an interest from 
the physical point of view for the following reason. While 
some of the qubit units are complicated devices, others, 
such as quantum dots, are produced at random. It might 
be the case that in the process of fabrication not one 
qubit unit was placed at some point, but a number of 
them very close to each other. Then, if there is no pos- 
sibility to directly study the structure of the chain for 
some reason, we provide simple tools to gain information 
about the system’s topography. 

Assume now, that we have a pseudo-chain with xa in- 
teraction between nearest blocks of spins, and local mag- 
netic fields acting on whole blocks. We start with the 
application of the Hamiltonian tomography by Burgarth 
and Maruyama ag, who showed that working in OES 
we can learn all the parameters of a linear chain Hamil- 
tonian could be learned from operating on the first site 
only. The observed dynamics would be as if the Hamil- 
tonian was 


N 
1 
Hy = 5 5 (Ji( Xi Xen + YiYi+1) + Bi — Z:))- (1) 


i=1 


Here with i = 1,2,...,N we number spin blocks (cur- 
rently there is only one spin in each block), (X;, Y;, Z;) 
is the vector of Pauli matrices acting on spin in the site 
i (except for i = N +1, where these operators vanish) 
and J;’s and B;’s are coupling constants, and local mag- 
nitudes of the magnetic field, respectively. The magnetic 
field terms have been chosen in a way that H|00...) = 0 . 
We will refer to the system described by as the model 
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FIG. 1: Above: a graph representing a pseudo-chain with 
three spins in the second block and two in the third. Dotted 
and dashed lines represent intra-block couplings K2 and K3. 
Below: the respective model chain with the same end-to-end 
dynamics in OES. 


chain. We assume that the actual system is described by 
N Ni Ni+i 

1 J; 
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x (XijXi+i,k + Yi jYi+1,k) 


Ni 
+ Bi (1 — Zij) 
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Here, N; is a number of spins-4 in the ith block, Ny = 
Ny = 1 and K; represents the in-block coupling. This 
interaction causes a modification of the effective magnetic 
field observed in OES, Bi = B; + (Ni — 1) Ki. 

One possible way to distinguish a pseudo-chain from a 
chain basing only on the at-ends data would be to inject 
two excitations, one from each end. The initial state is 
hence |y} = |10...01). Now, in a linear chain the two 
excitations propagate, might come close to each other, 
but in a pseudo-chain, they can come to the same block 
with more than one spin and pass each other, or then go 
back to the ends of the chains. The lowest number of 
applications of H on |y} to observe this effect is 2N — 
2. N — 1 of them allow excitations to reach the block 
and the other half lets them go back to the end spins, 
where they can be observed by the users. If the passing 
-1 
the excitations reach the rendez-vous point and as many 
to return. If there is N; spin in the block, there are 
again N;(N; — 1) distinct ways for the rendez-vous to 
happen-first choose to which spin goes the first excitation 
that reaches the block, then the second. The relative 
amplitude that the two excitation propagate along the 
chain meeting at this block is NP []j=1 J7- Additionally, 
either of the two excitation could be the first to leave the 


happens at block 7, there are ways in which 


meeting point, which produces extra 2 factor. Combining 
all these factors we reach the conclusion that 
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Of course, as energy is a global property, it cannot be 
learned directly from observing two extreme spins. How- 
ever, this difference has its effects on the system evolu- 
tion. Observers should deduce modified dynamics from 
the probability that after time t the system returns to 
state |w), 


P(|) + E) = || exp(LHt) |b)? (4) 


(throughout the paper we take h = 1 and J? = —1). Let 
us now expand the unitary evolution to the Taylor sum 
with respect to time. If we compare it to the evolution 
of |Y) governed by Ho, we obtain that 


IPU > Tele) = Pale) + WI = GG 
N-1 ; ~{N;-1 
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Our first procedure to reduce the class of possible geome- 
tries is as follows. First, one of the observers in turns 
brings the chain to the fully magnetized state, injects 
an excitation to find it back after some time, performing 
the Hamiltonian tomography (1g). It is not necessary, 
but certainly beneficial if two observers do it to average 
the results. Basing on their data, the users can now sim- 
ulate the behavior of two excitations in a linear chain, 
and measure their actual propagation. If they are able 
to identify the difference of the order of t?’~?, they could 
fit the possible nets. In simple words, one can construct 
a table of theoretical values of the right-hand side of 
in function of N;’s, and choose these, which are closest 
to in modulo, but below the value deduced from the ex- 
periment. The underestimation of the first term in the 
Taylor sum follows from the fact that higher order terms 
in the sequence will contribute with the opposite sign, 
as the probability function has an oscillative character. 
Note that the procedure described above does not allow 
us to learn intra-block coupling constants K;. Neither 
is it possible to order values of the pair {N;, NN—i+1}- 
These questions may be, at least partially, resolved in the 
other tomography routine. Burgarth and Maruyama 
showed how to conclude coupling constants and effective 
local magnetic fields starting from a pure state, and this 
remains the first step of our procedure. Nevertheless, it 
was proved by Di Franco, Paternostro, and Kim uA, that 


for a linear xx chain, this is equally possible with any ini- 
tial state, in particular, the maximally mixed state. This 
is obviously not true for a pseudo-chain, as this feature 
relies on the linearity of the Hamiltonian in the fermion 
picture. Hence we can use this difference to determine 
some features of the pseudo-chain. 

In a linear chain, such as our model (I), the initial 
operator has the following form 


X(t) = e 
= a (t)Xı + a2(t)Z1Y2 
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+ Q3 (t)ZiZ2X3 +... 
+ 61(t)¥1 + b2(t)Z1X2 
+ 3 (t)Z 122Y3 +... 
(6) 
with a; = Tam FTÌLI bi = = Da 0 Sig and 
bij = (—1)* (Ji-16i-1,5-1 + Jiði+1,j—1 
+ Biyij-1); (7) 
ing = CD (Ji1i-z1,j-1 + Iiiz, j-1 
E Bi: 68) 


The initial conditions read ð; o = 1 if i = 1 and 0 other- 
wise, and yio = 0. 

Assume now, that the first block containing more than 
one spin is labeled as i. Let us take (1 + +X) 
as the initial state. The non-trivial part of the state 
reaches the block after at least i — 1 commutations with 
the Hamiltonian. We now focus on the evolution of 

[]j<:_1 Z4) Xi-1, where we take i even; when it is odd, 


Xi-—1 is replaced with Y;_1, but the line of the argument 
remains the same. Apart from propagating back to the 
beginning of the chain, but might also propagate through 
the hub to spins in the block. This part of evolution leads 
to 
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j<i-l 


Ni 
XO Yar. (9) 
k=1 


First let us consider the effect of the magnetic field. In 
the model, where we have observed the effective magnetic 
field B; we shall have 


Ni 
(haa Z) D kei Yi 
> =1(Bi + (Ni) Ks) (jce-a 25) ERs Xar + (10) 
which in next i — 1 steps of evolution can retract to the 
beginning of the chain and contribute to (X1(0)Yi(t)) = 
{1(t). When we study the actual chain, however, only the 


real magnetic field B; allows the operator to contribute 
to the auto-correlation function. Hence the we have 
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where the auto-correlation function is normalized to 1 at 
t=0. 

Now, let us consider the evolution composed of 2i 
steps. There are a few trajectories such, that after the 
total of 2i commutations with the Hamiltonian the oper- 
ator comes back to the form of X1. The first is that the 
operator propagates back. This has a perfect analogue 
in the linear chain. Then it may propagate to (i + 1)th 
site/block and return following a different path than the 
first part of propagation, which obviously impossible in 
the model chain. Also, the operator propagated to the 
block can be twice a subject to the magnetic field and the 
intra-block exchange. The operator might also retract to 
the (¢ — 1)th spin and again propagate through one of 
N; branches of the hub. This is again in full agreement 
with the model evolution. The extra possibilities arise 
from the fact that e.g. Zi—-1X;,ı does not commute with 
Xj-1Xi,2 + Yi-1Yi,2. The propagatee gets a component, 
which is “double-headed” and proportional to 


(Taa Z) 
Ni Ni 
x (Yin era Xi kXi l + Xi-1 eA Yin Xi) 


Then either of the two “heads” can be the first to retract 
to the origin. After analyzing all these trajectories we 
come to 
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In the bracket the first term corresponds to the scenarios 
with the double-headed operator, the second arises from 
the operator passing the block but going back the other 
way, and the rest originates from the difference between 
the real and the effective magnetic field in the block. 

The only two quantities unknown from first routine are 
N; and K;. At most a finite set of possible solutions to 
Eqns. and can be found, especially if we assume 
that the coupling strengths are positive. 

The second procedure to estimate the structure of the 
pseudo-chain is the following. Assuming that we know 


the effective coupling constants J; and magnetic fields Bi 
we simulate the dynamics governed by the model Hamil- 
tonian (i). Next, we bring the system to maximally 
mixed state. Then we measure X on the first qubit ob- 
taining result a = +1. After some time tm,» we measure 
Xı getting bs = +1 in each run. We repeat the mixed 
state initialization and the measurement to reconstruct a 
few points of the mean function, fe(tm, s) = (abz(tm,x))- 
We then subtract these values from those calculated for 
the model chain, find a polynomial in time fit to the 
difference and put it to Eqn. (£2). Then a similar sub- 
routine is done, but with Y; as the second measurement 
and the result is used in the Eqn. (I). It is possible 
to measure Xı + Yı, but then the precision of the fit is 
compromised, as the function must be decomposed into 
odd and even parts. After solving the pair of equations 
we might update our model by introducing the first block 
to it, simulate the new dynamics and by comparing with 
the actual data recognize the structure of next block. Fi- 
nally, one can use the two-excitation procedure to verify 
N;,’s. 

One should mention that because the system cannot 
be mapped to a free fermionic field, the state before the 
Xı measurement cannot be arbitrary, but must be max- 
imally mixed. One must eliminate the possibility of hav- 
ing an operator initially acting on spins in a block, which 
would then evolve to X: or Yj. 

A question arises if the knowledge about the defects 
can help us in the task of the state transfer. Consider 
a pseudo-chain bathed in an environment, which causes 
local bit flips. We assume that the coherence time is 
much longer than the transfer time (whether in the state- 
mirroring chain , [l6] or in the bucket scheme (2H). 
This allows the high fidelity of the transfer, but at the 
same time causes a difficulty, when dealing with a pseudo- 
chain. When the initial state was the fully magnetized 
one, with probability A the excitation generated by de- 
coherence in the ith block will be able leave it and reach 
to one the ends, where it could by taken out by one of the 
users (when they constantly upload the |0) state). But 
in fraction & i= of all cases it is residual to the block. 
In this sense, blocks act as traps. It changes the mag- 
nitude of a spin by which the block could be effectively 
replaced and a coupling to it. The dynamics in the sin- 
gle excitation space becomes now an incoherent mixture 
of two evolutions, characteristic for two different linear 
chains: the one described by Ho and the one effectively 
modified by the trapped excitation. One can now calcu- 
late such an instant of time after which the probability 
for the excitation injected at the beginning to be trans- 
ferred to the last spin is 0, but it is sharply between 0 
and 1 for the model chain. If we measure Zy and get +1 
(corresponding to |0)), we exclude neither possibility. If, 
however, the measurement reveals an excitation transfer 
at this time, it could have been only due to the evolution 
without the excitation trapped in block 7. Successively, 


we can “empty” more traps in this fashion, one at the 
time. 


In summary, we present procedures, which allow to lo- 
calize additional spins in a non-linear pseudo-chain with 
zx interaction. Importantly, our method requires only 
data available from measurements performed on the ex- 
treme qubits and the control over a global parameter, 
which we associate with temperature: at its value 0 it 
should bring the system to the maximally magnetized 
state, at infinity it introduces the maximal entropy. Not 
only do we make it possible to estimate the number of 
spins in each block of the system, but we are also able to 
tell the intra-block coupling strengths. 


We stress that our procedures reveal different features 
of the system than those described in {14}. Therein, the 
system is still studied in OES, which allows it to have 
“hidden” states, i.e. a subspace with trapped excitations. 
It is, however, a very complex problem to to combine our 
procedures with ones of f 

It remains an open question whether there are com- 
putation algorithms, which can be performed more effi- 
ciently with pseudo-chains than with chains. Certainly, 
while xx chains are useless for computations, the dis- 
cussed systems might find interesting applications in this 
direction. The reason is precisely this: non-linearity 
of the structure introduces interactions betweens quasi- 
particles. We leave these possibilities for future investi- 
gation. 

Finally, we argue have argued that once we possess the 
knowledge about the structure of the pseudo-chains, we 
can remotely eliminate residual errors, which can appear 
due to locally acting decoherence. This allows to achieve 
almost perfect fidelity thanks to procedures such as the 
bucket scheme id. It is not possible without the removal, 
when the evolution is effectively a statistical mixture of 
evolutions of two different chains. 
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